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Summary 

The  primary  purpose  of  this  contract  was  to  assist  the  Naval  Ocean 
System  Center  in  implementing  a  fault  diagnosis  algorithm  previously  developed 
under  ONR  contract  by  the  principle  investigator.  To  this  end  our  major 
activity  was  directed  towards  the  development  of  an  efficient  data  base  for 
the  symbolic  transfer  function  data  required  by  the  algorithm.  The  resultant 
data  base  is  described  in  detail  in  this  report  and  is  characterized  by  both 
reduced  storage  requirements  and  reduced  retreival  time  as  compared  to 
previous  approaches.  By  combining  this  data  base  with  the  ATE  and  interface 
work  carred  out  by  NOSC  we  believe  that  a  viable  fault  diagnosis  algorithm 
for  linear  analog  circuits  can  be  implemented.  The  remaining  step  is  to 
actually  encode  a  software  package  which  incorporated  these  ideas. 

I .  Introduction 

Historically,  symbolic  network  analysis  has  been  motivated  by  the  prob¬ 
lems  of  circuit  design  and,  as  such,  the  emphasis  has  been  placed  on 

quickly  and  efficiently  obtaining  symbolic  transfer  function  from  a 
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given  set  of  circuit  specifications.  ’  In  an  operational  or  maintenance 


environment,  however,  one  is  typically  given  a  prescribed  nominal  circuit 
and  desires  determine  the  effect  of  various  (possibly  large)  perturbations 


thereon.  This  is  the  case  in  a  power  system  where  one  is  given  a  fixed 
network  and  desires  to  determine  the  effect  of  proposed  modifications  thereto. 
Alternatively,  in  the  problem  of  analog  circuit  fault  diagnosis  one  desires 
to  simulate  the  effect  of  a  number  of  alternative  failures  to  compare  the 
simulated  data  with  the  observed  failure  data.^ 

In  such  an  operational  or  maintenance  environment  numerous  perturbations 
of  the  nominal  circuit  are  studied  and,  as  such,  significant  computational 
efficiencies  can  be  obtained  if  one  first  generates  a  data  base  in  terms 
of  the  nominal  circuit  parameters  and  then  extracts  the  appropriate  symbolic 
transfer  function  from  the  data  base  each  time  a  different  symbolic  transfer 
is  required.  Of  course  the  benefit  to  be  achieved  via  such  an  approach  is 
dependent  on  the  size  of  the  data  base  and  the  ease  with  which  a  symbolic 
transfer  function  may  be  retrieved  therefrom. 

The  obvious  manner  in  which  to  generate  such  a  data  base  is  to  simply 
pre-compute  the  coefficients  of  all  required  symbolic  transfer  functions 
and  store  them  in  the  data  base.  Retreival  from  such  a  data  base  is,  of 
course,  immediate  but  the  data  base  may  become  overly  large.  Indeed,  the 
number  of  transfer  functions  which  must  be  stored  is  O(k^)  where  k_  is  the 
total  number  of  potentially  variable  circuit  parameters  and  £  is  the  maximum 
number  of  circuit  parameters  which  may  vary  simultaneously.  An  alternative 
approach  is  to  store  the  nominal  transfer  function  information  and  then  use 
Householder's  formula^  to  compute  the  required  symbolic  transfer  functions. 

In  such  a  data  base  we  need  only  store  0(£  )  transfer  functions  where  £  is 
the  total  number  of  component  output  terminals  but  retreival  requires 
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0(n  +p  )  multipl ications  where  p  is  the  actual  number  of  circuit  parameters 
which  vary  simultaneously.  Since,  in  practice,  ri  >>  p  the  retreival  process 
requires  approximately  0{n  )  multiplications  and  is  dominated  by  the  large 
dimensional  matrix  multiplication  required  by  Householder's  formula  rather 
than  the  low  dimensional  inverse. 

In  the  present  report  we  will  formulate  an  alternative  data  base  for 

o 

the  symbolic  transfer  functions  which  also  requires  0(rf )  entries,  but  for 
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which  retreival  requires  only  0{p  )  multiplications.  Since  p  is  typically 
small  this  is  tantamount  to  immediate  retreival. 

In  the  remainder  of  this  introduction  we  will  review  the  properties  of 
the  component  connection  model  for  a  large  scale  circuit  or  system^  which 
serves  as  the  starting  point  for  our  theory.  The  data  base  and  retreival 
formulae  for  the  case  where  p  _<  2  are  formulated  in  section  2.  while  the 
general  retreival  formula  is  derived  in  section  3.  Section  4.  is  devoted 
to  the  problem  of  retreiving  sensitivity  formulae  from  the  data  base  while 
section  5.  deals  with  the  problem  of  updating  the  data  base  when  the  nominal 
circuit  parameters  are  changed.  Finally,  section  6.  is  devoted  to  an  example 
illustrating  the  theory. 

The  component  connection  model  is  an  algebraic  model  for  an  inter- 
connected  dynamical  system  which  subsumes  the  classical  topological  models 
but  is  more  readily  manipulated  both  analytically  and  computationally.  The 
motivation  and  justification  of  the  model  are  discussed  in  detail  in 
reference  1  and  will  not  be  repeated  here.  The  component  connection  model 
takes  the  form  of  the  set  of  simultaneous  equations 

b  =  Z  ( j  lo)  a  1.1 
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1.2 


a  =  L^b  +  L^u 


y  =  L^-jb  +  L22u  1  -3 

Here,  Z  (=Z( ju>) )  is  a  frequency  dependent  matrix  characterizing  the  decoupled 

system  components  with  composite  component  input  and  output  vectors  a  and  b, 

respectively.  On  the  other  hand  the  L..;  i,j  =  1,2;  matrices  are  frequency 

i  j 

independent  connection  matrices  characterizing  the  coupling  between  the 
composite  component  vectors,  a  and  b,  and  the  Composite  system  input  and 
output  vectors,  u  and  y,  respectively. 

A  little  algebra  with  the  component  connection  equations  will  readily 


reveal  that 


S  =  L22  +  L21(1  -  ZL^r'zL 


IV  12 


where  S  (=  S(jw))  is  the  composite  system  transfer  function  matrix1 
characterizing  the  external  behavior  of  the  system  via 

y  =  S( juju  i_5 

Often,  rather  than  working  with  the  entire  S  matrix  we  find  it  convenient 

to  work  with  its  individual  entries;  sqv,  q  =  1,2 . £  and  v  =  1,2,  ...,v: 

which  are  related  to  the  component  connection  model  via 

sqv  =  Lqv  +  Lq  M  71  v 

L22  L2V‘  "  ZLlV  ZL12  1.6 

Here  Lq*  is  the  q-v  entry  in  l_22;  q  =  1,2,  ...  ,£  and  v  =  1,2,  . . . ,  v ;  is 

the  qth  row  of  L 2y  q  =  1,2,  ...  ,£;  and  L^2  is  the  vth  column  of  L12; 
v  =  1 ,2,  ...  ,  v. 

Finally,  since  we  are  interested  in  analyzing  the  effects  of  perturbing 
one  or  more  components  from  their  nominal  values,  we  decompose  Z  into 


nominal  and  perturbation  terms  in  the  form 


where 


Z  -  Zo  +  Z1 


|  cVrk 
k=l 


1.7 


1.8 


k  k  k  k 

Here,  c  (=  c  (jw))  is  a  column  vector,  r  {=  r  (jw))  is  a  row  vector,  and 

It 

6  is  the  scalar  perturbation  for  the  kth  potentially  variable  component 

k  k  k 

parameter.  In  a  typical  application  one  is  given  c  ,  r  ,  and  6  ; 

k  =  1,2,  ...  , k^;  characterizing  J<  potentially  variable  component  parameters 
though  at  most  £  such  parameters  vary  in  any  given  analysis;  p  k_. 

Indeed,  p  <<  J<  in  most  applications.  Finally,  we  note  that  can  be  expressed 
more  concisely  in  the  form 

Z  =  CAR  1.9 


where 


cP3 


1.10 


and 


1.11 


The  above  described  notation  formulated  for  the  component  connection  model 
is  summarized  in  table  1. 
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t 


Matrix 


Type 


Dimension 


1  ndex 


a 

composite  component  input 
vector 

m  x  1 

b 

composite  component  output 
vector 

ri  x  1 

u 

composite  system  input  vector 

v^  x  1 

y 

composite  system  output  vector 

a  x  1 

Ln 

connection  matrix 

m  x  n 

l21 

connection  matrix 

£  x  n 

iq 

l21 

qth  row  of  L21 

1_  X  £ 

q  =  1 

L12 

connection  matrix 

m  x  v^ 

L12 

vtb  column  of 

m  x  1 

v  =  1 

L22 

connection  matrix 

a  x  v 

$ 

q-v  entry  in 

1  x  1 

q  =1,2, 

s 

composite  system  transfer 
function  matrix 

a  x  x 

sqv 

q-v  entry  in  S 

1  X  1 

q=l,2,  .. 

z 

composite  component  transfer 

n  x  m 

Zo 

nominal  composite  component 

n  x  m 

Z1 

composite  component  transfer 
function  perturbation  matrix 

in  x  m 

ck 

column  vector  characterizing 
perturbation  of  kth  parameter 

n  x  1 

k  =  1, 

c 

k 

array  of  the  c  vectors  for  the 
parameters  chich  actual  vary 
(  row[ck3) 

n  x  p 

rk 

row  vector  characterizing 

1  x  m 

k  =  1, 

R 

l< 

array  of  r  vectors  for  the 
parameters  which  actually 
vary  (col[rk3) 

p  x  m 

6k 

kth  variable  parameter 

1  x  1 

k  =  1, 

A 

k' 

array  of  6  's  for  parameters 
which  actually  vary  (diag[6k3) 

P  x  p 

Table  1.  Summary  of  Component  Connection  Model 
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• •  •  » Q 

*  2 ,  •  •  •  ,  v 

*  •  •  »£>  v— 1 , 2  >  ...  i v 

*  -  »a»  v— 1  *2,  ...  i v 

2  k 

2,  ...  ,k 


2,  ...  ,k 


II.  The  Data  Base 


Our  data  base  is  composed  of  the  following  family  of  (frequency 
dependent)  scalar  transfer  functions 

So  =  L22  +  L21^  '  ZoLll>  ZoL21  ;  q  =  1  »2»  **•  v  =  1,2 . v  2.1 

b  L21^  "  ZoLll)  cJ  ;  q  =  1,2,  ...  ,£;  j  =  1,2,  ...  ,k_  2.2 

dkv  -  rkn  ♦  L„n  -  VnC’v-u'  k  ■  >-2-  -  •*  *  ■  ’-2-  -  *  2-3 

and 

e  J  =  r  Ln(l  -  ZqL^)  cj  ;  k.j  =  1,2,  . . .  ,k  2.4 

Here,  £  and  £  denote  the  number  of  external  system  inputs  and  outputs 
which  are  typically  few  in  number.  As  such,  the  ek^  array  composed  of  k_Z 
entries  dominates  the  data  base.  Also  note  that  all  of  the  entries  in  the 
data  base  are  formulated  in  terms  of  the  nominal  component  values  and,  as 
such,  the  data  base  may  be  generated  off-line  without  a  priori  knowledge 
of  the  perturbations  to  be  analyzed.  Finally,  the  entire  data  base  may  be 
generated  with  the  aid  of  only  a  single  £  by  £  (sparse)  matrix  inverse. 

Now,  if  we  assume  that  only  a  single  parameter  is  perturbed,  i.e. 

Z1  =  ck6krk  2.5 

for  some  fixed  k  =  1,2,  ...  ,k^,  to  retrieve  sqv  from  the  data  base  we  must 
evaluate 

sqv  =  Lq*  +  -  CZ0  +  ck6krk]L11)"1[Zo  +  ck6krk]L^  2.6 

*  k 

in  terms  of  the  elements  of  our  data  base  and  the  variable  parameter,  6  . 

To  this  end  we  invoke  Householder's  formula^ 
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2.7 


(W  +  XY)'1  =  W"1  -  W_1X(1  +  YW",X)'IYW 
with  W  =  (1  -  ZqL -j ^ ) *  X  =  -ck6k,  and  Y  =  rkLj.j  obtaining 
(1  -  [Zo  4  cWDL^r1  =  [(I  -  Z  L  )  4  (-ck6k)(rkL  )]_1 


-lvx-lw,rl 


•  0-Z0L„)-'  *  (l.Z0l1|)-|C,l'l|l-4,Lnll-;oL]|)-|c,(V'rlL]1|HoL]]) 


r1ck6ko-  -k- 

4-1  k.k  k, 


u  oir  ;  xm 


1-  IS  F 

-  6  e 


2.8 


Now,  upon  substitution  of  2.8  into  2.6  we  obtain 


cQV  _  .  qv  ,  .q  /,  ._k  k  k  \  “  1  r7  j.,.k  k  k  V 

S  “  L22  +  4l'  o  c  6  r  ■'hr  Cz0+c  6  r  •1Ll 


12 


'  L22  +  L21(,-ZoLlirlcZ»+cksl<rk:ILV 


12 


L21(1-Zohl)"lck6krkLii(1-2-Lii)'1tz^+ckfikrk]Li 


irl  o  i i 


12 


1  -  6kekk 

A^qk  k.  I  1  7  I  \~l-7  lV  j.  /rk^.qk  kk  k.  V 
.  -QV  .  .k.qkjc.v  .  6  b  r  hlf'-Zoh,)  Z0L12  +  (O  f  q  r  L 

-  so  +  6  b  r  Li2 +  - : — nrn - 

1  -  6  e 


■c 


skbqkdk,  +  (sk)2t.bqkeJ<krkLV2  +  ^k^v^  ^ 

i  i  ■  —  S  ^  ■  i 


,  .k  kk 
1  -  6  e 


so  *  j  ~  xmr z-9 

i  -  o  e 


which  is  the  desired  symbolic  transfer  function. 

If  we  assume  that  two  parameters  are  perturbed;  that  is 


Z,  =  cW  +  cW  2.10 

a  similar  formula  can  be  obtained  wherein  Householder's  formula  is  applied 
twice.  Since  this  formula  is  subsumed  by  the  general  retreival  formula  derived 
in  the  following  section,  we  simply  state  the  result  without  proof.  In 
particular, 
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III.  Retreival  Theorem 


% 


As  is  apparent  from  equation  2.11,  our  retreival  formulas  are  quite 
complex,  even  for  the  case  p  =  2  and,  as  such,  a  more  compact  notation  is 


required  if  they  are  to  be  tractable.  To  this  end,  we  assume  that 


•6  ;  k  =  1,2,  ...  ,p:  denote  the  potentially  variable  parameters  and  that 


Z.  =  l  ck<5krk  =  C  A  R 
1  k=l 


Of  course,  the  same  expression  applies  to  any  set  of  p  potentially  variable 
parameters  given  an  appropriate  change  of  the  index  set.  To  obtain  the  re¬ 
quired  symbolic  transfer  function  for 


s  =  42+L2l!’-tW-n>  t2o+Zl]L12 


with  the  above  specified  Zj  we  now  define  the  following  matrices  made  up  of 
elements  from  our  data  base 
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12 

$lv 

0 

so  ••• 

0 

21 

s22  ... 

s2^- 

qi  «.q2 

sSv 

’o  0 

0 

11  .12 
>  b  ... 

b]p 

21  .22 
i  b  ... 

b2p 

a1  t-a2  ... 

bflp 

V 


H 


|_dpl 

dP2 

dp^J 

4 

'e11 

e12  ... 

elp“ 

E  = 

21 

e 

22 

e 

e2p 

LeP’ 

ep2  ... 

epp_ 

while  A  is  defined  as  per  equation  1.12. 

THEOREM:  Using  the  above  notation 

s  =  L22+L21(1-[Z0Z1]L11)-1[Z0+Z13L12  =  So  +  B(  1  -  ae)']ac 


Proof:  First,  we  observe  that 

So  =  L22  +  L21(1_ZoLll^  ZoL12 
is  just  the  nominal  system  transfer  function  matrix  while 

6  ’  L2i<1-zoLnr'c 


and 


D  =  R[  1  +  l,1«1-Z0L1i)'',Z0]L1z  -  R{1-I.11Z0)-1 


via  Householder's  formula.  Finally, 


E  =  RLI1(1-Z0L11)-,C 


where  R  and  C  are  as  defined  by  equations  1.10  and  1.11.  As  such. 


(1  -  aE)'1  =  (i  -  arl110-zol11)"1c)-1 

-  [1  +  ARLn(l  -  (l-ZoL11)"1CARL11)'1(l-ZoLn)'1C] 

=  [1  +  ARL^l  -  ZoLll  -  Z1Lnr1C]-_ 

=  [1  +  ARL11(1-ZL11)“1C]  3.10 

where  we  have  invoked  Householder's  formula  with  Z  =  1,  X  =  aRL^,  and 
Y  =  (1-ZU.j  )~^C;  and  equation  1.9.  As  such, 

so  +  b(i  -  ae)_1ad  =  s0  +  l21{i-zcjl11)"1cci  +  arl11(i-zl11)"1c]ar(i-l11zo)_1l12 

*  VL2i<>-zoLn>-'tzi  *  zi4i<1-ZLn>'V<1-Liizo>% 

■  so  +  L2i<1-zoLn)''(t('-ZLii>*ziLii1('-ZLn)'1,2i(1-Ln2o>'1h2 

■  sotL2i<,-zoLii>'1<’-zoLnH'-2L,i>'lz)('-Luzo)'lLi2  - 

'  So  +  L2i(1-zl1i)"'Z1<1-|-11Z0)'1l12 
‘  L22  +  L21(1'ZoLlir'ZoL12  +  4l(,-ZLll>'lzi<1-LnZo)'lL12 

"  L22  +  L2lV‘L1lV  L12  +  L21*1'ZL11*  lzl*’'LnZo*  L12 

=  L22  +  L21fZot  "-ZLll)'lzi](,-LllZo)''ll2 
-  L22  +  L2i(1'ZLn)1[(,'ZLn)zo  + 

•  l-22  +  4l<,-ZLll>"1tZ  -  ZL11Zo]<'-L11Z1>''L12 
"  L22  +  L2l(1-ZLll)',zn-LnZo3n-L11Z0)-1L12 

“  L22  +  L21^'ZLn*  'ZL12  "  S 


as  required.  HU  Hill 
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IV.  Sensitivity  Formulae 

If  one  is  working  directly  with  the  component  connection  model,  it  is 
well  known ^  that  the  sensitivity  of  S  with  respect  to  a  parameter,  61,  can  be 
computed  via  the  formula 

■L2i<1-ZLn>'1[fi]n  +  Ln(1-ZLn)'1z:ILi2  ,1-1 

and  hence  it  is  appropriate' to  ask  whether  or  not  such  a  sensitivity  matrix 
can  also  be  computed  from  our  data  base.  Since  the  expression 

S  =  SQ  +  B(1  -  &E)_1aD  4.2 

is  formally  identical  to  1.4,  if  1  <_  i  <_  p  we  may  write 

+  E(1-4E)‘'o]D  4.3 

Ld«’J 

where 


dS 

dfi1 


M. 

l 


1 


0 


4.4 


with  the  one  appearing  in  the  ith  diagonal  entry.  Clearly,  the  expression 
can  be  computed  directly  from  the  data  base  with  the  same  level  of  computational 
effort  as  required  for  the  retreival  formula. 

In  the  case  where  61  is  not  included  in  the  given  set  of  parameters  which 
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deviate  from  nominal,  i  >  p  in  our  notation,  we  must  first  augment  the  B,  E,  D, 
and  A  matrices  to  include  6^  and  then  apply  equation  4.3  to  the  augmented 
system.  To  this  end  we  let 


B1  = 


D1  = 


E1  = 


b11 

b12  ... 

blp 

b11 

b21 

b22  ... 

b2P 

b2i 

& 

t£2 

b 9P 

b*1 

i11 

d12  ... 

d^~ 

i21 

d22  ... 

d2^ 

,pl 

d^2  ... 

dP^- 

11 

di2  ... 

d^_ 

11 

12  ... 
e 

Ip 

e 

e11 

21 

e22  ... 

2p 

e  K 

e2i 

Pi 

e^2  ... 

epP 

epl 

il 

elZ  ... 

e1-P 

ii 

e 

4.5 


4.6 


4.7 


and 


A3  - 


'  0 


4  . 

_ I _ _ 

o  i  o 


4.8 


The  we  obtain  the  retreival  formulae 


S  =  S  +  B1(l  -  AaE1)'1AaD1 


and 


ras 

dfi1 


=  B1  (l-AaEi)'1M  i.[l+Ei(l-AaEi)"1AaDDi 


p+1 
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4.9 


4.10 


V.  Updating  the  Data  Base 

In  many  applications  one  uses  a  data  base  such  as  that  described  above, 
as  a  design  tool  to  aid  in  simulating  the  effects  of  various  proposed 
modifications  to  the  system.  When  such  a  modification  is  finally  implemented 
it  is  then  necessary  to  update  the  data  base  to  reflect  the  new  nominal 
parameter  values 

Zo  ’  lc  4  J,  ‘W  •  z0  4  C4R  5-’ 

with  the  aid  of  Householder’s  formula  we  may  compute 

<'-zoLn>  -  tO-ZoM,)  -  wl,,]-1  - 

4  ('-z0Ln)-'cn  -  aRL11(i-zol11)-1c]-1arl11(i-z0l11)-1 

'  (1"2oLn)"1  4  5.2 

which  upon  substitution  into  equation  2.4  yields 


r  kl  k2 
[e  e 


ekp](l-AE)"1A 


Similarly, 

bqJ  =  bqk  +  [bql  bqZ  ...  bqp](l-AE)"1A 


5.3 


5.4 


dkv  =  dkv  +  [ekl  ek2 


ekp](l-AE)_1A 


5.5 
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and 


iqV  =  sqv 
0  0 


+  [bq1  bq2  ...  bqp](l-AE)-1A 


As  such,  the  entries  in  our  data  base  can  be  updated  with  a  computational 
effort  which  is  commensurate  with  that  required  by  the  retreival  formula. 


VI.  Examples 


Consider  the  simple  RC  op-amp  circuit  shown  in  figure  1.  The  component 


connection 


Figure  1:  RC  Op-amp  circuit. 


model  for  this  circuit  takes  the  form 


sC  0 


OR  0  i 


0  0  K  v 


0  -1  -1  i 


1  vi 


*  1  0  0  +  0 
'T 

0  1  0  Vo  1 


■•]■!!  "  ]  '•!  ‘  ["]  ■<] 


Thus  if  all  components  taken  to  have  nominal  values  of  1  we  obtain 


1  s  s 

(1"Z  J-n  )  =-110 


0 


63  0  0  1 


6.11 


c363r3 


0 

1 


k  i 

Combining  the  appropriate  c  and  r*  matrices  with  the  above  expressions  as 
per  equations  2.1  thru  2.4  we  obtain  the  data  base 


srt  =  1  6.12 

o 


b1  =-s 

b2.-i 

b3  =  s+1 

6.13 

d1  =  1 

d2  -  0 

d3  •  1 

6.14 

and 

en=0 

e’2  =  0 

e13  - 

e21  -s 

e22  =  0 

a23  -  -s 

e31  -  -s 

e32  -  -1 

33 

eJJ  =  s 

6.15 

where  we  have  deleted  the  q  and  v 

indices 

since  we  are  dealing  with  a 

single- 

input  single-output  system. 

— 

Now,  if  one  desires  to  compute  the  symbolic  transfer  function  with  respect 
to  perturbations  in  the  op-amp  gain  we  have 


,  .3  3  .3 

S(s.«3)  - 

0  1-6  e 


1  +  6_ 

1  -  63s 


6.16 


3 

Recalling  that  6  represents  a  perturbation  from  a  nominal  parameter  value  of 

3  3 

Kq  =  1  our  actual  gain  is  K  =  Kq  +6  =  1+6  ,  which  upon  substitution  into  6.16 

yields 


S(s,K) 


K 

(l-K)s  +  1 


19 


6.17 


which  is  the  classical  gain  formula  for  such  a  circuit. 

Finally,  if  we  desire  to  update  our  data  base  to  reflect  a  new  nominal 
value  for  the  circuit  parameters  of  C  =  1 ,  R  =  1 ,  and  K  =  2  we  invoke 

3 

equations  5.3  thru  5.6  with  6  =  1  yielding 


.3.3,3 

c  4.  _ ^ _ -  = 

So  ,  .3  33  . 

1  ~  6  e  63=1 
n  J3,3  31 

_11  ,  e  6  e  _  ^  _l 

e  + : — X33  ‘ 0 
1  -  6  e 


6.1 

6.1 


and  similarly  for  the  other  elements  of  the  data  base. 


VII.  Conclusions 

The  preceeding  development  has  been  motivated  by  operational  and  main¬ 
tenance  considerations  rather  than  the  design  considerations.  In  such  an 
environment  one  typically  deals  with  a  fixed  nominal  system,  but  carries 
out  repeated  analyses  thereon.  As  such,  the  cost  of  generating  the  required 
data  base  is  secondary  compared  to  the  cost  of  storing  the  data  base  and 
retreiving  information  therefrom.  In  these  respects  we  believe  that  our 
data  base  is  near  optimal.  Since  the  number  of  system  inputs  and  outputs 

p 

is  typically  small  our  data  base  contains  approximately  Ik  elements 
2 

(actually  k^  +  kjy+£)  +  vc^)  where  is  the  total  number  of  parameters 
which  are  potentially  variable.  This  data  base,  however,  contains  sufficient 
information  to  permit  one  to  retrjftve  symbolic  transfer  functions  for  any 
number  p  <_  k  of  variable  parameters.  Indeed,  the  number  of  variable  para¬ 
meters  in  a  symbolic  transfer  function  is  reflected  only  in  the  cost  of 
retreival  which  is  on  the  order  of  p  multiplications  (actually  p  +  p  _v  + 
pv(cj+l)).  Since  p  is  typically  small,  say  five  or  less,  this  is  minimal. 
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